Superfluid Weight vs. Superconducting Temperature based on a U(l) 
Slave-Boson Approach to the t-J Hamiltonian 
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Abstract 



Based on an improved U(l) slave-boson approach to the t-J Hamiltonian, we investigate 
a relationship between the superfluid weight n s /m*(the superconducting charge carrier 
density/the effective mass of the charge carrier) and the superconducting temperature 
q 1 T c . From the present study we find a linear increase of n s /m* with T c with the doping 



concentration in the underdoped region, a saturation around the optimal doping and a 
decrease in both n s /m* and T c in the overdoped region. Such a trend of the 'boomerang' 
shaped locus in n s /m* vs. T c with increasing doping concentration from the underdoped 



in 

£^ ■ to the heavily overdoped region is predicted to be in complete agreement with muon- 

o . 

spin-relaxation measurements. The boomerang behavior is found to occur in correlation 
with reduction in the spin singlet pairing(spinon pairing) order in the heavily overdoped 
region. 
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The transverse field muon-spin-relaxation(/i-SR) measurements of the magnetic penetration depth 
A in high T c copper oxide superconductors reveal an universal linear increase of the superfluid weight 
n s /m* (the superconducting charge carrier density / the effective mass of the charge carrier) with the su- 
perconducting transition temperature T c in the underdoped region, a saturation near the optimal doping 
and a decrease in both n s /m* and T c as the hole doping concentration increases from the underdoped 
region to the heavily overdoped region [JJ || 0. In the underdoped region the observed linear increase 
in both n s /m* and T c scales well with the doped hole concentration p while n s /m* vs. T c does not scale 
with p near and over the optimal doping. One can conjecture that the observed reduction of n s /m* 
with increasing hole carrier concentration in the overdoped region is attributed to either a decrease 
in the superconducting charge carrier density n s or to an increase in the effective mass of the charge 
carrier m*. In the present study we report a study of the observed boomerang behavior || [[| in n s /m* 
vs. T c as a function of doped hole concentration by using an improved U(l) slave-boson theory which 
correctly allows coupling between the spin and charge degrees of freedom. Contrary to our earlier U(l) 
slave-boson approach M , such introduction of both the spin and charge degrees of freedom was found 
to predict the arch shaped bose condensation temperature as a function of doped hole concentration p 

0, in agreement with the experimentally observed phase diagram in the plane of T c and p. For this 
study we first compute the doping and temperature dependences of the superfluid weight and discuss the 
cause of reduction in both the superfluid weight and the superconducting temperature in the overdoped 
region. In addition, a numerical comparison with the SU(2) slave-boson theory is briefly made to discuss 
a difference. 

We write the t-J Hamiltonian, 

H = -t ( c l c j* + cc.) + J ^ (Si ■ Sj - -UiTij) -fM ^2 i CiCT ' ( X ) 

<i,j> <i,j> i,a 

where Sj • Sj — \n(nj = — |(c|i c ]'T — c ll c ]l)( c j1 Ci l ~ c j! c *t)- Here S; is the electron spin operator at site 

1, Si = lc\ a cr a pCij3 with cr a p, the Pauli spin matrix element and rii, the electron number operator at site 
i, rii = c\ a Ci a . Using the single occupancy constraint and thus Cj CT = b\f i a (with f ia , spinon annihilation 
operator of electron spin a and b\, holon creation operator at site i), the U(l) slave-boson representation 
of the above t-J Hamiltonian leads to Eq.(Q), 

H = -t ( flfj^i + cc.) 

<i,j> 
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<i,j> 



-Mo 



E /U* + 2 E + 6 & - *)> ( 2 ) 



where Aj is the Lagrange multiplier field to enforce the single occupancy constraint at each site. In 
the SU(2) slave-boson theory, the electron operator is given by c a = -j=bSil) a with a = 1,2, where 

M ( h \ M 

tpi = and ^2 — an d h = are respectively the doublets of spinon and holon 

annihilation operators in the SU(2) theory. The SU(2) slave-boson representation of the above t-J 
Hamiltonian shows 



H = -\ E 

<i,j>a 

+(4/.i)(4&ii-4^2i) 

+ (/2i/l i -/u/2i)(fe t li^ + &Lfe2i) 

+(/iV 2 t l -4/i t ,)(4A, + 4/H) 

J 
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■^(l-ft)(l-^)x 



(A/L- - flfh){hjhi - Mu) - mo E ^ 

- E + + iAP^/m/m + blbu) 

i 

+^? ) (flifu ~ hifl + blbu - 4«] , (3) 

where ^M 2 ^ 3 ) are the real Lagrangian multipliers to enforce the local single occupancy constraint in 
the SU(2) slave-boson representation ||. It is noted that in both the U(l) and SU(2) approaches above, 
coupling between spin and charge degrees of freedom is correctly introduced in the Heisenberg term. 

After Hubbard Stratonovich transformations for the direct, exchange and pairing channels, we find 
the total free energy in the functional integral representation, for the U(l) slave-boson theory ||, 

F(A) = --In I D x DA f DA b DbDf 
P J 

e -(s b (A, x A f A b ,b)+sf( x A f J)) ^ ( 4 ) 
where S b represents the holon action for the charge degree of freedom, 
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S\A, X ,Af,A b ,b) = ^ dr\j2bH^r)(d T -^ b )b(r t ,r) 

Jo i 

+^Ei4i 2 K-i 2 +. 2 ] 

<i,j> 

-t J2 e iAij Xijtf( r iiT)b(r j: T) + c.c. 



<i,j> 



J 



(5) 



<i,j> 



and is the spinon action for the spin degrees of freedom, 



S'(*,A',/) = /%T[^/t(r„r)(a T -ii')Uh,t) 



<ij> 



J 



t(! -P) 2 J] XijfK r u T )f<r(rj,r) + c.c. 



<i,j> 



J 



(1-P) 2 E A£(/ i (r i ,r)/ T (r i ,r)-/ T (r i ,r)/ i (r i ,r))+c. C . 



(6) 



<«J> 



Here b is the holon field and /, the spinon field. Xiji A{j and ■ are the Hubbard Stratonovich fields 
corresponding to the exchange, the spinon pairing and the holon pairing channels respectively. After 
integration over the holon and spinon fields, we obtain the total free energy, 



F(A) = ~^J D X DAfDA b e-^ b ^ Af A b )+Ff( xA f)) } 



(7) 



where F b (A, x, A f , A b ) = -± In / £)6 e -s i (A )X ,A / ,A»,6) is the holon free energy an d F^( X ,A f ) = 
— ^ In J Dfe~ sf( - x ' Af 'f\ the spinon free energy. 

The linear response of current to weak applied electromagnetic(EM) field is, in the energy-momentum 
space, 



jl(u,q) = -Ui m (u,q)A m (u,q), 



(8) 



with j, the current and A, the EM vector potential. Here the current response function n ;m (u;,q) is 
obtained from 



Il lm (u},q) = -(3 



d 2 



dAi(-u, -q)dA m (u } q) 



F(A) 



(9) 



A=0 
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with (3 = k^p, where F(A) is the free energy of Eq.(^) or Eq.(??). The superfluid weight ^ is defined 
as the transverse EM current response function in the static, long wave length limit |§, 

Ti I 

^ = -limU xx (cu = 0,q), (10) 

for the isotropic system. 

As is shown in the second term of Eq.(^), the EM vector potential field A is seen to modulate(is 
coupled with) the hopping order parameter field Xij m the U(l) slave-boson representation. The phase 
fluctuations of the hopping order parameter, i.e. Xij = X e%aij are introduced. For the hole pairing of 
d-wave symmetry we allow the d-wave spinon pairing order parameter, Aj { = ± Ay (the sign +(— ) is for 
the ij link parallel to x (y)) and the s-wave holon pairing order parameter, A^ = A&. The EM current 



response function of the total system is obtained from the Ioffe-Larkin composition rule ITT 



IPV,q)IP>,q) 

%,q =^ r, 11 

n b (w,q) + n/(^,q) 

where n 6 (c<;,q) is the holon current response function to the the EM field A and IF (u, q), the spinon 
current response function to the gauge field a corresponding to the phase fluctuations of hopping order 
parameter. The holon and spinon response functions are calculated from the linear response theory(see 
appendices 1,2). The holon response function is evaluated to be 

nL(w,q) = [-2*xcosMAqi) 2 "Wqi)) - AqO V(-£ 6 ( qi )))<W 

qi 

+4t 2 X V^ sin(g lm + Q f) sin( qu + |) x 

{y (q + qi) v (qi) 2 - u \ q + cay (cay (q + (ch)) gg^^ g^g 

+(-« 6 (q + qi) V( qi ) 2 + M b (q + qi)^(qi)^(q + qi)^(qi)) 



n b (E b (q + qi)) - n b (-E b (q 1 )) 
^-(^(q + qO+E^qx)) 

+( _ Aq + + Aq + qi)tAqiK(q + <h y (qj)) i^+^^) 

+(Aq + qi) V(<h) . - .>(, + «a)u we + 



:i2i 



where -E 6 (q) = y (e fc (q) — A^) 2 — (</A 2 A?,7 q ) 2 is the holon quasiparticle energy and e 6 (q) = — 2tx7 q , the 
single holon quasiparticle energy with 7 q = (cosg^ + cosg y ). n b (E) = efJ E_ 1 is the boson distribution 



function, u b (q) = -^J e + 1 and v b (q) — 6 j^qj^ 1- The spinon response function is 
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qi 



+ (^(1 - Vfxf^ sin(g lm + *±) sm(q u + f ) x 

//p// , ^ ^(q + qi)(^ + ^(q + gi)) + (e f (g + gi) - n f )(e f (<a) - tf) + A /(q + qi)A^(gi) 

. f( f( Ef( qi )(-iuj + g^qQ) + (e% + qi) - fJ. f )(e f (<li) - V f ) + A^(q + qQA^qQ 

^ 1 J E/(q 1 )[(-^ + J E;/(q 1 )) 2 - J E;/(q + q 1 ) 2 ] 

_ f f Ef(g + qiK-jo; + g^q + qQ) + (e% + qj - //)(e%i) - //) + A^(g + gi)A^( qi ) 

[U 1 lq + qijj J E/(q + q 1 )[(-^ + E/(q + q 1 )) 2 -E/(q 1 ) 2 ] 



, u f( ^(qi)(^ + ^(qi)) + (6% + qi ) - ^)(e/( qi ) - //) + A^(q + qi)A^( qi ; 
^ 1 E/(q 1 )[^ + E/(q 1 )) 2 -E/(q + q 1 ) 2 ] 



(13) 

where E^(q) = J (e^ (q) — /i-^) 2 + (Aj(q)) 2 is the spinon quasiparticle energy, e^(q) = — f (1 — p) 2 X7q; 
the single spinon quasiparticle energy, and A'^(q) = J(l — p) 2 Af(p n , the spinon pairing gap with <£> q = 
(cosg x — cosg y ). n*(E) = e $E +1 is the fermion distribution function. For the U(l) theory the superfluid 
weight is obtained from Eq.(|T0D with the use of (|TT|). In the present calculations of the superfluid weight 
J/t = 0.2 is chosen. For other choice of J/t we find that there exist no qualitative differences in the 
behavior of n s /m* vs. T c . 

In Fig.l the computed superfluid weight displays a linear increase in doped hole concentration(rate) 
p in the underdoped region(the predicted optimal doping rate is p Q = 0.07). As hole doping concen- 
tration further increases, the superfluid weight first saturates and rapidly drops to in the heavily 
overdoped region. This behavior in the overdoped region is attributed to a decrease in holon pairing 
interaction caused by the predicted diminishing trend of the spinon pairing probability amplitude (order 
parameter) A^-, as can be readily understood from the last term of Eq.(||). The predicted decrease 
of pseudogap(spin gap) temperature with hole doping rate is in excellent agreement with observation. 
In Fig. 2, a "boomerang" shaped locus in " s< -^°' > and T c is displayed by showing a linear relationship 
between the two in the underdoped region and the 'reflex' behavior in the overdoped region. Although 
not numerically agreeable, this trend is completely consistent with muon-spin-relaxation measurements 
[0]. As mentioned above, the spinon pairing amplitude A^- is further reduced in the overdoped region 
and this, in turn, causes a decrease in holon pairing interaction and thus the superconducting charge 
carrier density n s . This will allow the reflex behavior, that is, a decrease in both the superfluid weight 
and superconducting temperature. 

For the sake of introducing the low energy phase fluctuations of the order parameters Xij an d A^, 



which were not considered in the above U(l) slave-boson study, we now introduce the SU(2) slave-boson 
theory. In the SU(2) theory, the EM field A is found to modulate both the hopping order parameter Xij 
and the spinon pairing order parameter A{-. The phase fluctuations of both Xij and A^- are thus taken 
into account to assess the EM current response of the system. We introduce the spinon pairing order 
parameter of d-wave symmetry, Ajj = ±A/(the sign +(— ) is for the ij link parallel to x (y)) and the 
holon pairing order parameter of s-wave symmetry, A^. a/3 = A b (5 ajl 5g ; i — S^Sp^)- 

In the SU(2) slave-boson theory, the response function II^, A ' a -* of holon isospin current to the EM 
field A and the gauge field a vanishes owing to the contribution of the 6 2 -boson in the static and long- 
wavelength limit (see the appendices 3,4). Therefore the superfluid weight of the total system is given 
by the holon current response function only, 

^ = in?rV = o, q = o). (14) 

Here iT^f'^cj, q) is computed from the use of the usual linear response theory for the holon action 
Eq.(??)(see Appendix 6). In Fig. 3, we show the doping dependence of the superfluid weight with the 
choice of J/t = 0.2. The superfluid weight predicted from the SU(2) theory increases faster than the 
U(l) case, as doped hole concentration. This is because the gauge fields(or the phase fluctuations of 
the order parameters) do not screen the EM field in the SU(2) theory. In Fig.4, the relation between 
T c and ns{ \^^ is displayed; at low doping, T c increases linearly with ns ^T°^ as in the U(l) theory, a 
plateau near the optimal doping p = 0.13) and a reflex in the overdoped region is also observed. There 
exists only a quantitative (but not qualitative) difference between the two theories. This is because the 
spinon pairing amplitude A^- predicted by the SU(2) theory slowly diminishes in the overdoped region 
compared to the U(l) case. 

In summary, we investigated a relationship between the superfluid weight and the superconducting 
temperature, that is, vs. T c based on both the U(l) and SU(2) slave-boson approaches to the t-J 
Hamiltonian. From both theories we find a qualitatively similar boomerang shape in the path of " s 
vs. T c by showing a linear increase in the underdoped region and a reflex behavior in the overdoped 
region. Such trend of boomerang behavior in high T c superconductors is completely consistent with 
/i-SR experiments. The reflex behavior predicted by both theories is attributed to the diminishing 
attractive hole pairing interaction caused by the markedly reduced spin pairing order A^, particularly 
in the heavily overdoped region, and thus to the reduction of both the superconducting charge carrier 
density n s and the superconducting temperature T c in the overdoped region. 
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FIGURE CAPTIONS 



Fig. 1 The doping dependence of the superfluid weight in the U(l) slave-boson theory at T/t = 
0.001 with J/t — 0.2. p is the hole doping rate. 

Fig. 2 The plot of Ws( - T ^ ; 001f - ) vs T c in the U(l) slave-boson theory with J/t = 0.2. The thick arrows 
point towards increasing doping rate. p Q = 0.07 is the predicted optimal doping rate. 

Fig. 3 The doping dependence of the superfluid weight ^ in the SU(2) slave-boson theory at T/t = 
0.001 with J/t = 0.2. p is the hole doping rate. 

Fig. 4 The plot of " s(r =°; 001f) vs T c in the SU(2) slave-boson theory with J/t = 0.2. The thick arrows 
point towards increasing doping rate. p Q = 0.13 is the predicted optimal doping rate. 



9 



FIGURES 
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